Depinning of a dislocation: the influence of long-range interactions 
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The theory of the depinning transition of elastic manifolds in random media provides a framework 
for the statistical dynamics of dislocation systems at yield. We consider the case of a single flexible 
dislocation gliding through a random stress field generated by a distribution of immobile dislocations 
threading through its glide plane. The immobile dislocations are arranged in a "restrictedly random" 
. manner and provide an effective stress field whose statistical properties can be calculated explicitly. 

' We write an equation of motion for the dislocation and compute the associated depinning force, 

which may be identified with the yield stress. Numerical simulations of a discretized version of 
the equation confirm these results and allow us to investigate the critical dynamics of the pinning- 
' depinning transition. 
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I. INTRODUCTION 

The depinning transition of a dislocation moving in a glide plane has been investigated in the past in several 
papers usinga line tension approximation for the dislocation self-interaction and considering interactions with point 
obstacles jlB. In fact, however, the bending of a dislocation line produces a long-range self-stress [HQ and the 
obstacles responsible for the pinning are often forest dislocations rather than point defects (i.e. solute atoms, small 
clusters). This distinction is important since dislocations, in addition to short-range interactions, give rise to long- 
q ' range stresses in the glide plane and thereby cause a long-range correlated pinning field. Also the 'short-range' 
interactions (junction formation) with forest dislocations which in real crystals contribute substantially to the pinning 
should not be modeled in terms of point obstacles, since the corresponding pinning forces exhibit spatial correlations 
with ranges of the order of the junction length, which is in turn of the same order as the obstacle (forest dislocation) 
■ spacing j|. 

Since long-range interactions are known to affect the general scaling properties of the depinning transition H , in this 
paper we study dislocation depinning considering explicitly the effects of nonlocal elasticity and spatial correlations in 
the pinning field. To this end, we introduce a lattice model with long-range dislocation self-interactions that correctly 
recovers previous results for the bending of a dislocation segment. Using this model we show that the geometrical 
properties of the dislocation at depinning differ from those expected in line tension models, as shown by a different 
roughness exponent. 
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We consider a dislocation gliding in the xy plane, parametrized by R(s) — (X(s), Y(s), 0), with Burgers vector 



g 

7 1 , II. FORMULATION OF THE MODEL 

c : 

, b = (b, 0, 0), and interacting with a pinning field created by a random array of immobile dislocations threading through 
the plane. The resolved shear stress acting locally on the dislocation is given by the sum of three contributions 

>< *UR) = + < Z {R) + <AR), (i) 

h ; 

5« | where a^f is the externally applied stress and a% z is the pinning field. We assume that electron and phonon drag lead 
to an overdamped motion of the dislocation, so that the velocity normal to the dislocation is given by v n = ba xz /T, 
where V is an effective viscosity coefficient. Using relations given by Hirth and Lothe (t]] for the stress fields of 
dislocation segments, the self-stress of the dislocation can be written as 

s ffo f (Y - Y(s))d s X + (X- X(s))d s Y/(l - v) 

where fi is the shear modulus and v is Poisson's number. 
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To study the influence of long-range correlations in the pinning stress, we consider a stress field created by an 
array of immobile dislocations threading through the glide plane. We note that it is not feasible to consider a fully 
random arrangment of the threading dislocations, as this would imply a diverging energy density of the pinning field 
H . Hence some correlations in the dislocation arrangement creating the pinning stress must be assumed to ensure a 
finite energy density. To this end we use Wilkens' construction of a 'restrictedly random' distribution of dislocations 
H where the system is divided in boxes of side I, each of them containing n dislocations with a fixed distribution of 
Burgers vectors placed at random within the box. For this model, it is possible to compute analytically the statistical 
properties of the pinning field. In particular, the stress distribution is symmetric with tails decaying as p/cr 3 , where 
p = n/l 2 is the density of dislocations creating the pinning field ||. The spatial correlation function of the pinning 
field can also be computed fl0|| . This function is given by 

KAR))°IAR')) = c{iibfpf{[R-R']/i s ) , (3) 

where the constant C depends on the distribution of Burgers vectors and £ s = M/^/p with M w 0.25-^/n is the 
correlation length of the pinning field. The correlation function /(i?/£ s ) has the following properties: / — ► — ln(R/£ s ) 
for R = \R\ — > and / = for -R/£ s > 4. The logarithmic singularity of / at the origin is truncated at the dislocation 
core radius r c ss b, yielding a mean square stress ((c^) 2 ) = p 2 b 2 Cp\n[^ s /b]. The spatial average of the pinning stress 
is, of course, zero. These results have been verified numerically. 

III. DEPINNING TRANSITION 

We consider a dislocation with average screw orientation (i.e., the average line direction is parallel to the x axis). In 
a first step, we perform a simplified analysis and assume that the dislocation has no overhangs (i.e R = (x,h(x),0)), 
and only small deformations (i.e. \h — h'\ <C \x — x'\). In this case, the equation of motion becomes 

T dh oxt bp f dx (iii 7/ <\\ (x — x') dh(x') 



bat ^ + fy^^^)-^0) + ^^j + ^M. (4) 

Eq. ^ describes an 'elastic' interface moving in a quenched disordered medium. For this class of systems we expect a 
depinning transition as a function of the external stress. The average velocity of the dislocation vanishes in the long 
time limit when the external stress a = c^f is below a critical value a c , and scales as (a — cr c )' 3 for a > a c . The 
critical point may be identified with the yield stress. 

We can obtain the dependence of the critical stress on the density of the dislocations creating the pinning field by 
rescaling Eq. (0)^ If we define a; = [4TTy/Cp\x, t = [AirC'ppb 2 /T]t and a = [a/ (pby/Cp\, the resulting equation depends 
only on [My/C] and the rescaled external stress, so that 



a c = b^y/Cp d c {u, [MVC]), (5) 

where a c is the depinning stress for the rescaled equation. When pinning is due to the cumulative effect of a large 
number of obstacles (weak pinning), this stress can be estimated using an argument which is originally due to Larkin 
On small scales and for weak pinningJJ], the self-interaction keeps the dislocation essentially straight. Hence the 
average force experienced by a segment of length L is y/ ((crP) 2 ) L£ s . Above a critical length L c the dislocation can 
bow and adjust to the 'contours' of the pinning field. Bowing a near-screw segment of length L over a characteristic 
distance £ s -C L requires a characteristic force ~ pb(l + z/)£ s /[(l — v)L\. Equating this to the pinning force, one 
obtains the 'Larkin length' L c = [pb[(l + v)/{l - v)} y/£, s / ({<r p ) 2 )} 2/3 ■ The depinning str ess can be envisaged as 
the strength of the pinning field averaged over the Larkin length, i.e. a c = y/ ((cr p ) 2 )£ s /L c . This estimate yields 
~a c tt[MVC{l-v)/(l+v)] l l\ 

Close to the depinning transition the dislocation is expected to be rough and satisfy a self-affinity transformation 
h(bx) ~ b^h(x), where £ is the roughness exponent. In practice, we can obtain £ from the power spectrum of h(x), 
which scales as S(k) ~ fc~( 2 ^ +1 '. Power counting suggests that the equation for the dislocation dynamics is in the 
same universality class as an equation in which the self stress is replaced by V 2 /i, which corresponds to a line tension 
approximation. Numerical simulations indicate that in the latter case the lineshape is superrough with £ ~ 1.25 [ [l2|| . 
In addition, transverse motion can lead to a fractal lineshape for strong disorder pj. 



1 In the context of dislocations weak pinning corresponds to the case of Mott-Labusch statistics Jll| . 
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IV. LATTICE MODEL 



In order to test the analytical considerations we use a lattice model which allows for transverse motion of dislocation 
segments in the x direction. We discretize the system into a grid of square cells with boundaries parallel to the x 
and y axes and we assign a "spin" Si = 1 to the slipped region and Si = — 1 to the unslipped region. The dislocation 
is defined as the interface between slipped and unslipped region. By definition, this implies that it is discretized in 
terms of screw and edge segments. 

The stress on a given dislocation segment is computed by summing up the external stress, the pinning stress and 
the stress due to the other segments which is given by 



where fi is the point of gravity of the ith segment, the first sum is over the screw segments and the second over the 
edge ones, <r s = S s pb 2 x/[4:Tr\r\ 3 ] and a e — S e pb 2 y /[4n(l — v)\r\ 3 are the stress fields due to screw and edge segments 
in the origin and the signs ^s, S e depend on the orientation of the segments. Periodic boundary conditions are 
imposed in the x direction by performing the Ewald sum. In order to avoid possible problems with the interaction of 
near segments , we identify the segment length with the dislocation core radius and neglect the self-energy of the 
segments. 

The dynamic evolution is performed by chosing at random one segment and moving it one step by flipping one 
spin in the direction of the local force. At each time step the stresses are updated. In order to test the effects of 
the random update and of the discretization of the dislocation, we considered the bowing of a dislocation line pinned 
between two points (Fig. 1). This configuration has been investigated in a numerical study by Foreman J3j. In spite of 
the severe simplifications we make, the critical stresses and critical configurations are in good agreement with results 
of that work. 

We simulate the lattice models using square lattices of size up to L = 1000. For the pinning field, we use again a 
random stress field created from a Wilkens construction performed by placing two screw dislocations of opposite sign 
and line direction normal to the glide plane at random in a box of size I. This corresponds to C — l/(47r), M = 0.35. 
In the simulation reported here we further assume v = 1/3 and measure stresses in units of p, and distances in units 
of b. For the dislocation density creating the pinning field, we use the value p = 0.00125. We obtain a c by studying 
the survival probability P s |lj] for small lattice sizes L ~ 100. In the limit L — > oo one expects P s {v) = 0(<r — <j c ), 
where 9(x) is the Heaviside function. For finite L, we impose that a c lies in the region for which 0.5 < P s < 0.7 and 
obtain a c = 0.0040 ± 0.0005 for p = 0.00125 (see Fig. 2). 

For comparison, we note that the Larkin estimate yields <r c = 0.0037. This value is consistent with the simulation 
result, although the Larkin length L c is not large as compared to the correlation length £ s as required by the "weak 
pinning" assumption. 

Simulating the model for larger system sizes, we see that the dislocation roughens in time but overhangs and islands 
are not observed (i.e. the lineshape is not fractal, see Fig. 3). We next compute the power spectrum of h(x) at <r c 
for a system of size L = 1024, averaging over ten realizations of the disorder, and find £ = 1 (see Fig. 4). This result 
differs from the value £ = 1.25 obtained using a line tension approximation |12|| . 



We have investigated the depinning of a dislocation moving in a fixed glide plane and interacting with a random 
field that has long-range correlations. For creating the pinning field, we have used a Wilkens construction with an 
array of immobile 'forest' dislocations threading normally through the glide plane. It is emphasized that this is not 
considered to be a realistic treatment of the interaction between moving and forest dislocations - we have neglected 
junction formation, while the Burgers vector and line direction of the pinning dislocations have been chosen in such 
a manner as to yield maximum stresses in the glide plane of the moving dislocation. Rather, it was our aim to create 
a pinning field which (i) has correlations with a range of about half the spacing of the pinning dislocations and (ii) a 
pinning strength which corresponds approximately to the strength of a random forest. (For the present parameters 
and using relations given in jl5j], this amounts to a c « 0.0060.) In this sense the present treatment is supposed to yield 
a more realistic picture of the depinning of a dislocation interacting with a dislocation forest than previous models 
assuming point obstacles |lj . 

As a result, we find that the use of a point obstacle field in conjunction with a line tension approximation tends to 
overestimate the roughness of the dislocation at the depinning transition. Long range self-interactions give rise to a 
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V. DISCUSSION AND CONCLUSIONS 
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roughness exponent which is smaller than the value obtained from a simple line tension approximation. The lineshape 
of a dislocation depinning from a dislocation forest is of significant importance for theories of work hardening and 
microstructure evolution. Fractal lineshapes of depinning dislocations may help to explain the observation of fractal 
dislocation patterns , and the fact that a moving fractal dislocation must leave behind loops has been invoked as 
a reason for work hardening Jl5| . While the present results do not support this viewpoint, further studies are required 
to analyze possible crossovers to fractal dislocation lineshapes in stronger pinning fields. 
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FIG. 1. Bowing of a dislocation pinned at the two ends, with at without side arms. For the two cases Ref. [1] reports 
g c ~ 0.0011 and a c ~ 0.0008 
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FIG. 2. The survival probability for L = 80, 120. 
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